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Abstract 

The density dependencies of various effective interaction strengths in the relativistic mean field 
are studied and carefully compared for nuclear matter and neutron stars. The influences of different 
density dependencies are presented and discussed on mean field potentials, saturation properties 
for nuclear matter, equations of state, maximum masses and corresponding radii for neutron stars. 
Though the interaction strengths and the potentials given by various interactions are quite different 
in nuclear matter, the differences of saturation properties are subtle, except for NL2 and TM2, 
which are mainly used for light nuclei, while the properties by various interactions for pure neutron 
matter are quite different. To get an equation of state for neutron matter without any ambiguity, 
it is necessary to constrain the effective interactions either by microscopic many-body calculations 
for the neutron matter data or the data of nuclei with extreme isospin. For neutron stars, the 
interaction with large interaction strengths give strong potentials and large Oppenheimer-Volkoff 
(OV) mass limits. The density-dependent interactions DD-MEl and TW-99 favor a large neutron 
population due to their weak p-meson field at high densities. The OV mass limits calculated from 
different equations of state are 2.02 ~ 2.8IM0, and the corresponding radii are 10.78 ~ 13.27 km. 
After the inclusion of the hyperons, the corresponding values become 1.52 ^ 2.06 Mq and 10.24 
~ 11.38 km. 

PACS numbers: 21.30.Fe, 21.60.Jz, 21.65.-hf, 26.60.-I-1 
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I. INTRODUCTION 



A widely used and successful approach for nuclear matter and finite nuclei is the mean 
field theory employing effective interactions. The mean field theory includes non-relativistic 
mean field theory with effective nucleon-nucleon interactions such as Skyrme or Gogny, and 
the relativistic mean field (RMF) theory. In a certain sense, the RMF theory is more funda- 
mental as it starts from a phenomenological hadronic field theory with strongly interacting 
baryons and mesons as degrees of freedom ill. It has been used not only for describing 
*e p.ope.*. of n„c.ei nea. the vaHey of JLy .ucce.sf„..y U. but a.o fo. predict., 
the properties of exotic nuclei with large neutron or proton excess Q] . In the mean field 
theory, the effective interactions are adjusted to various properties of nuclear matter and 
finite nuclei. In recent years, a number of effective interactions of meson-baryon couplings 
based on the RMF theory have been developed, including nonlinear self-couplings for the 



cT-meson and/or w-meson, such as NLl, NL2 ^, NL3 y, NLSH TMl and TM2 
However, these nonlinear interactions haveproblems of stability at high densities, as well as 
the question of their physical foundation jsf . A more natural alternative is to introduce the 
density dependence in the couplings P|. Based on the Dirac-Brueckner calculations, Typel 
and Wolter proposed the density-dependent effective interaction TW-99 and expected that 
the model could be reasonably extrapolated to extreme conditions of isospin and/or density 
j^. Along this line, Niksic et al. developed another effective interaction DD-MEl 10]. In 
this paper, we will analyze the density dependencies of various effective interactions includ- 
ing both the nonlinear and density-dependent versions in the RMF theory and investigate 
their influences on properties of nuclear matter and neutron stars. 

The existence of neutron stars was predicted following the discovery of neutron. In 1934, 
Baade and Zwicky suggested that neutron stars could be formed in "supernovae" The 
radio pulsars discovered by Bell and Hewith in 1967 jl^ were identifled as rotating neutron 
stars by Pacini ll3|l and Gold 14]. The flrst theoretical calculation of neutron stars was 

n n 

performed by Oppenheimer and Volkoff [l^, and independently by Tolman [la|. In their 
calculation the neutron stars were assumed as gravitationally bound states of neutron Fermi 
gas. Recent progress on the study of neutron stars can be found in Refs. jl^, I3| and 
references therein. 

The physics of neutron stars has offered an intriguing interplay between nuclear pro- 
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cesses and astrophysical observation, and has become a hot topic in nuclear physics and 
astrophysics. The neutron stars exhibit conditions far from those encountered on the earth. 
The neutron star models including various so-called realistic equations of state have resulted 
in the following general picture for the interior of neutron star. The surface of neutron star 
is a solid crust of thickness about 1 km, which is mainly made up of nuclei and free elec- 
trons. Inside the crust, charge-neutral neutron star mainly consists of neutrons together 
with a small concentration of protons and electrons in equal number. Proton and electron 
densities increase with total baryon density, and the vr, K mesons and other baryons 
( e.g., hyperons ) as well as a phase transition from baryon degrees of freedom to quark 



matter will appear 



The equation of state (EOS) is essential to understand the structure and properties 
of neutron stars. The EOS determines properties such as the mass range, the mass-radius 
relationship, the crust thickness, the cooling rate and even the energy released in a supernova 
explosion. Usually, the EOS is obtained by extrapolating the theory, which is developed 
mainly for normal nuclear matter, to nuclear matter with extreme high isospin and high 
densities. Unfortunately, such extrapolation is always model-dependent. The RMF theory 
has proved to be very successful in describing the properties of nuclear matter and finite 
nuclei jl|, Q], rotation nuclei nuclei far from (3 stability 0,0] and magnetic rotation jioj ]. 
etc. Based on the RMF theory, there are also many studies on neutron stars and strange 
nuclear matter with effective interactions including nonlinear self-couplings for scalar and 
vector mesons [17, 20-30]. 

In this paper, the density dependencies of various effective interaction strengths in the 
RMF theory are studied and compared for nuclear matter and neutron stars. The corre- 
sponding influences of different density dependencies for effective interactions are presented 
and discussed on mean field potentials, saturation properties for nuclear matter, EOS, max- 
imum masses and corresponding radii for neutron stars. In Section II, a brief description of 
the RMF theory in nuclear matter and neutron stars is presented. The results and discus- 
sions are given in the following section. In the last section, we give a brief summary. 
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II. A SKETCH OF THE RMF THEORY IN NUCLEAR MATTER AND NEU- 
TRON STARS 



The details of RMF theory can be found in a number of reviews |l|, |j, |17|. The RMF 
theory starts from an effective Lagrangian density with baryons, mesons {a, u and p) and 
photons as degrees of freedom [h = c = 1): 



B 

B 



1 — T"3^ 

h'^df, - TUB- gaBcr - g^Bl^^^^l - gpBl^TB ' - e'j'^A^ 



B 



+ - ^mla^ - U{a) - ^J^.uo^" + \mluo^uo>' + U{uo) 

where the Dirac spinor ^/^^ denotes the baryon B with mass vtlb and isospin r^. The sum 
on B is over protons, neutrons and hyperons (A, S^, S^, et. al.) in this paper. The 
scalar sigma (a) and vector omega [lo) offer medium-range attractive and short-range repul- 
sive interactions respectively, and the isospin vector rho (p) provides the necessary isospin 
asymmetry. Their masses are denoted by m^^m^ and vtip. The corresponding coupling 
meson-baryon constants are gaB-,9ujB and respectively, is the isospin of baryon B 
and is its 3-component. The nonlinear self-couplings for a and uj mesons are respectively: 

= ^g^a' + ^gsa\ U{uj) = ^C3(cuX)^ (2) 

with the self-coupling constants g2,g3, and C3. The field tensors u^^, p^^,, and A^^, are for 
w-meson, p-meson, and photon respectively. 

The parameterization of the interaction in the RMF theory is obtained by fitting the 
properties of nuclear matter and some finite nuclei. Instead of the self-coupling of the 
meson fields, Typel and Wolter proposed the density dependencies of the couplings in the 
RMF theory j^, i.e., the coupling constant ga{ui)B of the a {u) meson is replaced by: 

ga{u)B{p) = gaioj)Bipo)fa{c,){x), (3) 

where 

is a function of x = p/po) with the vector density p = \fj^i 3ii = 'n,B^Bl^l'^B and 
the saturation density of nuclear matter pg. The eight real parameters in Eq. ^ are not 



independent. The five constraints /cr(w)(l) = 1, f^(uj)(^) ~ ^ fai^) = f"i^) reduce 
the number of independent parameters to three. The density-dependent p-meson coupling 
constant gpB is introduced as, 

9pb{p) = 9pb{Po) exp[-ap(x - 1)], (5) 

with two parameters and gpB{po)- 

The equations of motion for baryons and mesons can be derived from the Lagrangian 
density in Eq. (^. In the following, we present only the case for the density-dependent 
couplings. The equations for nonlinear couplings can be obtained easily by adding the non- 
linear self-couplings of the mesons and neglecting the density dependencies of the coupling 
constants. 

The equations of motion for baryons are: 



i^B = 0, (6) 



where the effective mass m*^ = itlb + gaB<^ and is a "rearrangement" term due to the 
density dependencies of the couplings: 



Kb = '-fi^i^BY^BOO. + -f^i^B^TBi^BP, + -^V^bV^B^)- (7) 

The field equations for mesons and photons are respectively: 

{d'^dp + ml)a = -J2 g^rBi'Bi^B, (8) 

B 

dpu'^'' + mlu- = J2 Qu^Bi^BYi^B, (9) 

B 

d,p^^ + my = J2 apBi'BYTBi^B + gpBP, x P"", (10) 



B 



d.A^'' = e^^^^-^^. (11) 



Nuclear Matter 



For infinite nuclear matter, introducing the mean- field approximation, i.e., the meson 
fields are replaced by their mean values, and neglecting the coulomb field, the baryon wave 
function is the eigenstate of momentum A;, and the source currents '^b'^b and '^bI^'^b in 
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Eqs. fpiTUj) are independent of the spatial coordinate x. Thus the equations of motion can 
be simphfied as: 



[l^ikf, - g^B^^i, - QpeTb ■ Pt, - SJb) - ^*b\ i^B^k) = 0, (12) 

mla = -gaBps, (13) 

mluo = gujBp, (14) 

"^pP0,3 = ^gpBT3BpB- (15) 
B 



For the nonhnear self-couphng effective interactions, the corresponding terms — ?7'(cr) and 
—U'{uJo) should be taken into account in the Eqs. (fT^ and (fT^ respectively. 
The eigenvalues of the Dirac equation for baryons in Eq. (fT^ are obtained as: 

eB{k) = g^B^^o + QpB ■ r^B ■ Po,3 + ^qb + \/k'^ + m*j^, (16) 

where S^^ is the time component of the rearrangement term. The baryon vector density p 
and scalar density ps are respectively: 

B B B 

Ps = J:<^b^Pb>=J2psB = -,J2 I k'dk (18) 

B B B V +^B 

where fc^ denotes the Fermi momentum of baryon i?, and the no-sea approximation has 
been used. 

The energy density and pressure of nuclear matter are respectively, 

^ = \mW + + \mlph + k^dk^Jk^ + mf, (19) 

111 1 r'^^ k^ 

2 2 2 ^ ' ^ Stt^ ^ Jo V ^ + ^B 

B. Neutron Stars 

The charge-neutral neutron stars includes not only neutrons and protons, but also leptons 
A (mainly e~ and p~ ) in equal number to protons and also hyperons at high densities. The 



7 



Lagrangian density for neutron stars is similar to Eq. (^J), except an additional term for 
leptons: 

Yl M^Yd^-mx)^x, (21) 

X=e~ ,fi~ 

Introducing the mean-field and no-sea approximation, the equations of motion for baryons 
and mesons can be derived, and the corresponding energy eigenvalues, baryon density and 
scalar density can be obtained for neutron stars, similar to those for the nuclear matter. 
The equations of motion for electron and fi~ are free Dirac equations and their densities can 
be expressed in terms of their corresponding Fermi momenta as p\ = A;f/(37r^). 

The chemical potentials fiB for the baryons B are the energy eigenvalues of the Dirac 
equation: fis = ^b(^)- The chemical potentials fix for the leptons are the solutions of their 



equations of motion: nx = + m\. In neutron stars, the chemical equilibrium conditions 
are: 

/^B = bslin - QBlJ'e, IJ'fi = IJ-e, (22) 

where bs and Qb denote baryon charge and electronic charge of baryon B, /i„, fi^ and /x^ 
are the chemical potentials for neutron, electron and fi~ respectively. The baryon number 
conservation and charge-neutral conditions are given by: 

B B 

o = e<?b+5:«^=e#-e|^=o- (24) 

B X B X 

The energy density and pressure for neutron stars are respectively: 

1 22,-'- 22,-'- 22 



p = i:p^=i:i^' (23) 



1 



+ ^^y k^dk^k^ + m*j^ + — I k^dk^k^ + ml , (25) 

B ^ X=e~ ,fi~ 

D -'-22|-'-22|-'-22 , v^i 



^OB 



B 

1 V- f^^ k^ „ 1 — 



o V / dk + ^ V / dk , . (26) 
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III. RESULTS AND DISCUSSIONS 



For nuclear matter, Eqs. (fT ^ - (fTHj) provide a set of coupled transcendental relations 
defining the meson fields and energy eigenvalues. The list of unknowns is 

P, ks, Ps, QujB, gaB, gpB, 0", C^O, P0,3, £B{k). 

For a given baryon density p and the asymmetry of the nuclear matter t = {pn — Pp) / p, we 
can get the coupling constants g^B, fl'^s and gpB from Eqs. ©-(0), as well as the neutron 
density p„ and proton density pp, and of course their corresponding Fermi momenta fc„ and 
kp. The uj and p fields can be obtained by Eqs. ()14|) and ()15p and the a field can be solved 
from Eqs. (jl3p and (jl8|) by iteration. 

The properties of neutron stars can be obtained by solving the Eqs. flT ^ - (fTS|) and 
()24|1 by the following procedures: (1) For a given baryon density p, taking initial values of the 
meson fields (a, ujq, po.s) as well as the neutron and electron chemical potentials {pn, Pe), 
the particle densities and Fermi momenta for electrons {ke, Pe), protons {kp, pp) and 
hyperons {kh, ph) can be obtained via the eigenvalues in Eq. (fT^ . chemical equilibriums 
and charge- neutral conditions in Eqs. (j22j) and (j24|) . which in turn fix the neutron density 
Pn from baryon number conservations Eq. (2) With the particle densities (p^, Pe, and 

Pfj,) and Fermi momenta (fc^, ke, and A;^), the meson fields (a, ojq, p^s) can be obtained from 
Eqs. (fT ^ -(fT3 j) : (3) These two steps should be repeated by iteration until the self-consistence 
is achieved. 



A. Effective interaction strengths in nuclear matter and neutron stars 

Using the nonlinear RMF interactions NLl, NL2 Q, NL3 Q, NLSH Q, TMl, TM2 Q, 
and GL-97 17 1 and the density-dependent interactions TW-99 0], and DD-MEl 1^, the 
density dependencies of various effective interaction strengths in RMF theory are studied 
and carefully compared in nuclear matter and neutron stars. 

In Fig. n the density dependencies of the effective interaction strengths for a (top), 
UJ (middle) and p (bottom) mesons in symmetric nuclear matter as functions of the nu- 
cleon density are shown. The shadowed area corresponds to the empirical value of the 
saturation density in nuclear matter ( Fermi momentum kp = 1.35 ± 0.05 fm~^ or density 
p = 0.166 ± 0.018 fm~^). These curves are labelled from the top to the bottom at p=0.15 



fm~^ orderly from left to right. For the nonlinear effective interaction, the "equivalent" 
density dependencies of the effective interaction strengths for a, uj and p are extracted from 
the meson field equations Eqs. (fT^ - (fT^ according to: 

g<7B{,p) = 9aB + U'{a)/ps = gaB + + Q^iCT^) / ps, (27) 

gwB{.p) = Qi^B -U\ujq)/ p = g^B - {c^^^Dl (28) 
9pb{p) = 9pB- (29) 

The density dependencies of the interaction strengths for TW-99 and DD-MEl are very 
similar for symmetric nuclear matter in Fig. Q as noted in Ref. Here the comparison 
between the nonlinear and the density-dependent interaction will be emphasized. 

For the cx-meson, the interaction strengths given by TW-99 and DD-MEl are quite dif- 
ferent from the others in either magnitudes or slopes. In particular, strengths of TW-99 
and DD-MEl for the density interval in Fig. ^are almost twice as large as that of GL-97. 
Differences for nonlinear and density-dependent interactions can also be seen in the region 
of the empirical nuclear matter densities. For the w-meson, except TW-99, DD-MEl, TMl 
and TM2, all the other strengths are density- independent. All the strengths are similar 
to each other in the region of the empirical saturation densities compared with those of 
the (T-meson, although large differences can also be seen at low densities. For the p-meson 
which describes the isospin asymmetry, the strengths for TW-99 and DD-MEl show strong 
density dependencies in contrast with the constants in the other interactions. They cross 
the nonlinear interactions at a density much lower than the empirical saturation density. 

The interaction strengths as functions of baryon density for neutron stars matter are 
given in Fig. El As the effective interactions NL2 and TM2 are mainly used in light nuclei, 
we don't discuss them here. At densities p < 0.2 fm~^. Fig. |21is similar to Fig. ^ For the 
scalar a-meson, the interaction strengths of TW-99, DD-MEl, TMl and GL-97 decrease 
with the baryon density in similar slopes, while those of NLl, NL3 and NLSH decrease with 
baryon density for the densities p < 0.2 fm~^, then increase afterwards. This is due to the 
positive g^a^ in Eq. (j27p for NLl, NL3 and NLSH, in contrast with the negative ones for 
TMl and GL-97. For the vector cu- meson, the interaction strengths of TW-99, DD-MEl 
and TMl decrease with the baryon density. At densities 0.2 < p < 0.55 fm~^, the strengths 
of DD-MEl, TW-99 and TMl are between those of NL3 and GL-97. The curve given by 
TMl crosses with the line of GL-97 at density p ~ 0.57 fm~^, and then gives the weakest 
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interaction strength. For isospin-vector p meson, the interaction strengths of TW-99 and 
DD-MEl decrease with baryon density and trend to vanish at high densities, while the others 
are constants. 

Although the aspects of the interaction strengths are quite different from each other, as 
it will be shown in the following sections, all of them give fair descriptions for the properties 
of nuclear matter. 

B. Potentials for nuclear matter and neutron stars 

Potentials for symmetric nuclear matter, pure neutron matter and neutron stars calcu- 
lated with density-dependent interactions TW-99 and DD-MEl are illustrated in Figs. 01 
I3J and El The results are shown in comparison with those obtained with the nonlinear 
interactions, such as NLl, NL2, NL3, NLSH, TMl, TM2 and GL-97. 

As the contribution of p-meson potential for symmetric nuclear matter vanishes, we show 
the vector potentials q^jb^^o, scalar potentials gaB<y and their sum as functions of nucleon 
density p for different effective interactions which are marked from top to bottom at density 
p=0.15 fm^'^ orderly from left to right in Fig. jSl For the vector and scalar potentials, GL- 
97 gives the weakest results due to its weakest interaction strength in Fig. ^ The curves 
associated with the other interactions lie between those by TM2 and NL2. Their difference 
increases with the density, the difference of scalar potentials at p=0.15 fm~'^ between NL2 
and TM2 is around 120 MeV, and that of vector potentials is around 130 MeV. At saturation 
density, the scalar and vector potentials given by the density-dependent interactions TW-99 
and DD-MEl are similar to the nonlinear interactions except for GL-97, TM2 and NL2. 
For the total potentials for symmetric nuclear matter guiB^o+daBO', GL-97 also gives the 
weakest result and the difference for different interactions is about 20~30 MeV in the range 
of saturated densities. 

The vector potentials g^iB^o, scalar potentials gaBO', isospin potentials gpBPo,3 and their 
sum as functions of the nucleon density for different effective interactions are shown in Fig. 
mfor pure neutron matter. The variations of the vector potentials with density p are same 
as those for symmetric nuclear matter because they are related to the whole nucleon density 
p only. The variations of the scalar potentials with density p are slightly different from 
those for symmetric nuclear matter, as the same density for symmetric nuclear matter and 
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pure neutron matter implicates different Fermi momentum, fc^, as shown in Eq. psp . The 
p-meson provides the necessary isospin asjTumetry. In pure neutron matter, the isospin 
potentials gpBPo,3 given by nonlinear interactions increase linearly with the density due to 
their constant interaction strengths in Fig. ^ While for the same reasons, we can easily 
understand the isospin potentials for the density-dependent interactions. The compensation 
between the density dependencies of gpB and the increase of the density makes the isospin 
potentials increase at first and decrease after p = 0.12 fm~^. The sum of the vector potentials 
gojB^o, scalar potentials gaB'^, and isospin potentials gpBPo,3 gives the total potentials for 
baryon in pure neutron matter. The total potentials are attractive at low densities and 
become repulsive at high densities. Different from those of symmetric nuclear matter, due 
to the contribution of the isospin potentials for pure neutron matter, the difference between 
the total potentials is quite large, e.g., the largest difference 40 MeV between NL2 and 
TW-99 occurs at p = 0.15 fm'^. 

For neutron stars. Fig. El shows the vector potentials g^B^^o, scalar potentials gaBO', iso- 
vector potentials gpBPo,3 and the sum of them for different effective interactions as functions 
of the baryon density p. For scalar and vector potentials at densities p < 0.2 fm^^, the 
properties are similar to those in Figs. ElandlU The vector potentials g^B^o, which offer 
short-range repulsive interactions, increase with the baryon density, while the scalar poten- 
tials gaB<^i which offer medium-range attractive interactions, change evidently at densities 
p < 0.3 fm~^, and reach saturation at density p 0.4 fm~'^. Just as for nuclear matter, 
GL-97 always gives the weakest scalar and vector potentials. Due to the strong interaction 
strengths of NLl, NL3 and NLSH at high densities as shown in Fig. |2l they give large 
scalar and vector potentials at densities p > 0.3 fm"'^. The results calculated with DD-MEl, 
TW-99 and TMl lie in the middle. The difference in potentials between the two density- 
dependent effective interactions DD-MEl and TW-99 shows up clearly beyond the density 
p =0.4 fm"'^, in contrast with that for nuclear matter. 

As shown in Fig. El the contribution of the p-meson provides the necessary isospin asym- 
metry for neutron stars. For densities 0.065 < p <0.1 fm~^, the potential gpBPo,3 calculated 
with various interactions range from 20 MeV to 30 MeV. For nonlinear interactions, the 
potentials gpBPQ,3 increase with the baryon density and are about 100 MeV at density p=0.9 
fm~^, which are about 10% of the vector potentials g^jB^o for GL-97 and TMl, and 5% 
of those for NLl, NL3 and NLSH. For density-dependent interactions, the isospin poten- 
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tials trend to vanish after p=0.75 fm~^, due to the density dependencies of the p-meson 
interaction strengths in Fig. |21 

The total potentials QuB^o+QaBO' + gpBPo,3 are attractive at low densities. At high den- 
sities, they become repulsive and the neutron stars are bound by the gravity. The results 
given by various interactions are similar at densities p <0.2 fm^'^ and become quite different 
at p >0.3 fm~^. The total potentials for NLl, NL3 and NLSH are close to each other and 
increase rapidly. Those for TW-99, GL-97 and TMl are close to each other after p ~ 0.16 
fm~^, while that of DD-MEl lies in-between. 

C. Properties of nuclear matter 

In Fig. ini we display the energies per nucleon, Eb/A = e/p — m, for pure neutron 
matter (upper panel) and symmetric nuclear matter (lower panel) for density-dependent 
interactions TW-99 and DD-MEl as functions of nucleon density p. For comparison, the 
results for various nonlinear interactions are also shown. We label the curves from top to 
bottom at density p=0.20 fm~^ orderly from left to right. 

In the upper panel of Fig. IHl the energies per nucleon for all interactions are always 
positive and increase with the nucleon density for pure neutron matter. The results for 
density-dependent interactions DD-MEl and TW-99 are very similar to each other. At 
densities p <0.075 fm~^, they are larger than the other interactions except for NL2, thereafter 
they cross the curves for the other interactions and give the smallest energy per nucleon after 
densities p >0.17 fm~'^. 

For symmetric nuclear matter, the curves for energies per nucleon for various interac- 
tions display similar dependencies on densities below the saturation density and pronounced 
differences at higher densities. The saturation densities for various interactions are similar 
(p = 0.166 ± 0.018 fm-\ Eb/A = -16.33 ± 0.10 MeV) and are located in the empirical 
range of saturation density except those of TM2 and NL2, which are mainly used for light 
nuclei. 

From the upper and lower panels of Fig. IHl it can be seen that the equations of state at low 
densities (below saturation density) for pure neutron matter and symmetric nuclear matter 
are different. The deviations of energy per nucleon Eb/A for symmetric nuclear matter 
below saturation density are negligible compared with those for pure neutron matter (the 
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deviation is about 15 MeV at density p = 0.15 fm~^). The different results for symmetric 
nuclear matter and pure neutron matter mainly come from the potentials as shown in Figs. 
El and m especially from the isospin potentials gpBPo,3, which exhibits large deviation. 

From Fig. |21 we have seen that the difference for the total potentials at saturated density 
in symmetric nuclear matter is about 20~30 MeV, while the energies per nucleon at satu- 
ration density are close to each other in Fig. IHl To understand these, we will discuss energy 
per nucleon at the saturation density for symmetric nuclear matter in detail. 

From Eq. (fT^ and the meson field equations Eqs. (fT^ - (fT3j). the energy density for 
symmetric nuclear matter can be expressed as: 



1 r''^ 

B B 

-\gaB(yps - \g2<y^ - jg3(^^ - ^gc^BUJop + jc^u^ -^pb-^1 



^OB 



where 



+ [gujBt^OP + y2PB ■ ^OB + k'^dky/k'^ + {tub + gaB(yY ] 

B ^ B 

£a + euj+ £re + ^TV, (30) 



£a = ~gaB(rPs - ^g2(r^ " ^fi'sC^^, (31) 
1 1 4 

£uj = --gcuBUJop + -cst^O) (32) 

£re = -^Pfi'^OB' (33) 

B 

= g^B^Op + ^pB-^0B, (34) 

B 

^V/ k^dk^k^ + irriB + gaBaY, (35) 



k _ 
J — 



B 

£n = + ^AT) (36) 

in which, the contributions from a and uj fields are respectively and e^, the rearrangement 
term is Sre, and the contribution from nucleons is ^at with its potential energy part and 
kinetic energy part e\j. For nonlinear interactions, £:re=0, while for the density-dependent 
interactions, (72 = (73 = C3 = 0. 
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In TableHl the Fermi momenta kp (po = ^f/Svt^), vector densities po, scalar densities p^, 
scalar potentials (^o-bo", vector potentials Qujb^o, the various contributions €„, e^j, Sre, £%, 
and En, the system energy densities e, nucleon mass m and energies per nucleon Eb/A = 
e/p — m at saturation density for different interactions are shown. The saturation densities 
given by various interactions are similar (po ~0.150 fm~^), except for TM2 (po =0.132 fm"'^). 
From Eqs. (jHUI) - (jHUj) . we know that if the scalar densities ps and vector densities po (or Fermi 
momenta kp) at saturation densities are similar, larger scalar potentials (negative) will give 
larger Eo- (positive) and smaller e% (positive), and larger vector potentials will give larger 
(negative) and (positive). Taking into account the contribution of the rearrangement term 
Ere, the energies per nucleon Eb/A = e/p — m for various interactions become closer to each 
other. For example, for TW-99 and GL-97, which give the largest and smallest potentials 
respectively, the difference of QctbO' ~ -210.5 MeV leads to Asa =16.04 MeV fm~^ and 
Ae% =-30.80 MeV fm'^, and the difference of g^B^ ~ 194.1 MeV leads to Ae^ =-14.77 MeV 
fm~^ and Ae^ =29.70 MeV fm~^. Taking into account the contribution of the rearrangement 
term Ae^e =-0.16 MeV fm~^, the difference between the total energy densities is Ae =0.01 
MeV fm~^, i.e., the difference between their energy per nucleon AEb/A = A{e / pQ—m) =0.07 
MeV, which is small compared with Eb/A ^ 16 MeV. 

The saturation properties of symmetric nuclear matter for different interactions including 
the Fermi momenta kp, saturation densities po, energies per nucleon Eb/A, effective masses 
m* and m*/m, in compressibility K and symmetric energy coefficients agym are shown in 
Table |n] We can see that the results given by the density-dependent interactions DD-MEl 
and TW-99 are similar to those of the nonlinear interactions except for TM2, NL2 and 
GL-97. The saturation density for TM2 is lower and the energy per nucleon for NL2 is 
larger than the corresponding empirical values, which can also be seen in Fig. IHl Due to 
the weakest scalar potentials shown in the middle panel of Fig. 121 GL-97 gives the largest 
effective mass m*/m=0.78 and incompressibility K=2AQ MeV, which are justified from the 
empirical nuclear saturation properties 

Although properties below the saturation densities for symmetric nuclear matter are quite 
similar, the EOS at low densities given by various interactions for pure neutron matter is 
quite different. This is due to the effective interactions used so far are obtained by fitting 
the properties of doubly magic nuclei, which have an isospin close to that of the symmetric 
nuclear matter. To get an EOS for neutron matter without any ambiguity, it is necessary 
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to constrain the effective interactions eitlier by microscopic many-body calculations for the 
neutron matter data or the data of nuclei with extreme isospin. 



D. Properties of neutron stars 

The energies per baryon for neutron stars as functions of baryon density for different 
interactions are given in Fig. [7| At low densities, Fig. [7| is very similar to the top panel 
of Fig. El as the neutron stars matter is almost pure neutron matter at low densities. The 
energies per baryon for the nonlinear interactions NLl, NL3 and NLSH increase quickly with 
the density compared with those of TW-99 and GL-97, while the results for DD-MEl and 
TMl lie in-between. These results are in consistent with the total potentials for neutron 
stars in Fig. |S| At density p ~ 0.3 fm^'^, the potentials for TW-99 and DD-MEl are the 
lowest, while that of DD-MEl crosses GL-97 and TMl at p ~ 0.4 fm'^ in Fig. The Eb/A 
for TW-99 is the smallest. 

The evolutions of particle densities with the baryon density in neutron stars are given 
in Fig. ISJ and the corresponding figures in logarithm scale for low densities (0.05 < p 
fm^ < 0.2) are given in the sub-figures, where the solutions are from p=0.065 fm~'^ to 0.9 
fm-^ i.e. 0.425 < p^/po < 5.88 with po = 0.153 fm'^. 

At low densities, the charge-neutral neutron star matter is mainly composed of neutrons. 
As the density increases, high-momentum neutrons will /3-decay into protons and electrons ( 
n p + e~ +Vf, ) until the equihbrium at which the chemical potentials satisfy Pp = pn — Pe- 
As the neutron density increases, so do the proton and electron densities. When p^ attains 
the value of the muon mass, the p~ will appear. The equilibrium with respect to the reaction 

^ p~ +V^ + i'e implies that pe = p^- The p~ thresholds are different for different effective 
interactions. All the p~ thresholds are in the range of p = O.lliO.Ol fm~^ with the minimum 
and maximum thresholds given by NLl and GL-97 respectively. 

From Fig. |H1 it can also be seen that the density-dependent effective interactions TW- 
99 and DD-MEl give the largest neutron densities, and accordingly the smallest proton 
densities due to the baryon number conservation. This is because the strengths gpB for the 
density-dependent effective interaction become weaker with the baryon density, as shown in 
Fig. 121 Because of the charge-neutral condition, the densities of electron and p~ for different 
effective interactions have the same tendencies as the proton densities. 



16 



As hyperons would appear at roughly twice saturation density, it is necessary to sudy 



neutron star with hyperons. The details for t 
relativistic mean field theory are given in Refs. 



he mc 



One can introduce the ratios of the 



usion of hyperons in neutron star in 



meson- hyperons coupling strengths coupling constants to those of nucleons as: 

9(Th 9i^h 9ph ("^7^ 

9crN 9wN 9pN 

where, g^h, 9u}h and Qph have the same density dependencies as g^N, g^N and gpN respectively, 
and the ratios x^h = x^jh = Xph = a/2/3 are chosen according to Ref. 3^ . 

The equation of state (EOS) is very important to understand the structure of neutron star. 
The stiffer the EOS, the larger the mass that can be sustained against collapse. There are two 
constraints for the realistic EOS. One is a stiff limit by the causal constraint dp/ de < 1, which 
results in the limit mass of just over 3Mq. The other corresponds to the soft limit, which 
corresponds to the free Ferr ni g as with neutrons, protons and leptons in equilibrium and the 
limit mass is about 0.7 Mq |13]. Here the EOS calculated by different effective interactions 
are given in Fig. IHl The solid lines represent these with neutrons and protons only, and the 
dashed lines represent their corresponding ones with hj^erons included respectively. As can 
be understood from the potentials in FigEl the nonlinear interactions NLl, NL3 and NLSH 
give stiffer EOS than the other interactions, GL-97 and TMl give softer EOS, and density- 
dependent interactions TW-99 and DD-MEl lie in-between. Furthermore, the inclusion 
of hyperons softens the corresponding EOS considerably, as shown by the corresponding 
dashed lines in Fig. IHl The softest EOS is given by TMl. After the inclusion of hyperons, 
the corresponding solutions for the nonlinear interactions, NLl, NL3 and NLSH, exist only 
below density, p = 0.42, 0.51, and 0.58 fm~^, respectively. Beyond the corresponding density, 
the scalar potential gaB<^ will increase and make the effective masses m* = m+g„B<^ negative. 
Therefore in the following, we do not discuss the properties of neutron stars with hyperons 
for NLl, NL3 and NLSH. 

Fo. a .tat,c ,.oba> St., the OppenKei.e.Vo..off-To>„,a„ (OVT) e,„at,o. is QQ: 

dp _ [p{r) + e{r)][M{r) + ATTr^p{r)] , . 

dr ~ r(r - 2M(r)) ^ ^ 

M{r) =Att[ e{r)r^dr (39) 



The point R, at which the pressure vanishes, p{R) = 0, defines the radius of the star and 
M{R) is the gravitational mass. For a given EOS, the OVT equation has a unique solution 
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which depends on a single parameter characterizing the conditions of matter at the center. 
This can be chosen as the baryon density or energy density. In Fig. ^[ the masses versus the 
central densities (left panel) and radii (right panel) of neutron stars for different interactions 
are shown. The solid lines represent these with neutrons and protons only, and the dashed 
lines represent the corresponding ones with hyperons included. There is critical maximum 
value for the masses of neutron star, known as the OV mass limit. Beyond this mass the 
star is unstable to gravitational collapse. 

The OV mass limits, corresponding radii, central densities, energy densities and pressures 
for different effective interactions are presented in Table ITTTl The second rows for DD-MEl, 
TW-99, GL-97 and TMl represent the quantities with hyperons. Without hyperons, the OV 
mass limits calculated from different EOS are 2.02 ~ 2.8IM0, and the corresponding radii 
are 10.24 ~ 13.27 km. The OV mass limits and corresponding radii for DD-MEl and TW- 
99 are respectively 2A75Mq, 11.903 km and 2.195M0, 11.209 km and they give softer EOS 
and smaller OV mass limits than those of the effective interactions NLl (2.8O9M0, 13.137 
km), NL3 (2.778Mo, 13.081 km) and NLSH (2.8O3M0, 13.270 km). The results for TMl 
are similar to TW-99, while GL-97 gives the smallest OV mass limit (2.OI8M0) and radius 
(10.779 km). 

Due to the softer EOS for neutron star with hyperons, small OV mass limits have been 
obtained. The OV mass limits calculated with hyperons for DD-MEl, TW-99 GL-97 and 
TMl are 1.52 ~ 2.O6M0, and the corresponding radii are 10.24 ~ 11.38 km. With hyperons 
included, the OV mass limits and corresponding radii for DD-MEl and TW-99 are respec- 
tively 2.O6IM0, 11.375 km and I.868M0, 10.853 km. The results for TMl are similar to 
TW-99, and gives radius (11.366 km) and the smallest OV mass limit (I.517M0). 

IV. SUMMARY 

We have studied and carefully compared the density dependencies of various effective 
interaction strengths in symmetric nuclear matter, pure neutron matter and neutron stars. 
The corresponding influences on potentials and properties of symmetric nuclear matter, pure 
neutron matter and neutron stars are presented and discussed. As the interactions NL2 
and TM2 are aimed for light nuclei, we don't present their results for neutron star. The 
properties calculated by the interactions NLl, NL3 and NLSH are close to each other. The 



18 



same conclusion can be seen for the density-dependent interactions TW-99 and DD-MEl. 

For the cr-meson, all the interaction strengths are density-dependent. While for the u- 
meson, TW-99, DD-MEl and TMl are density-dependent. For the p-meson, the interaction 
strengths for TW-99 and DD-MEl decrease very fast with density while the others are 
constant. Even though the interaction strengths and the potentials from various interactions 
are different, the differences of saturation properties for various effective interactions in 
symmetric nuclear matter are subtle except for NL2 and TM2. 

Unlike those for symmetric nuclear matter, the properties for pure neutron matter by 
various interactions are quite different. As the effective interactions used so far are obtained 
by fitting the properties of doubly magic nuclei, it may be successful for nuclear matter 
with an isospin close to that of the symmetric nuclear matter. To get an EOS for neutron 
matter without any ambiguity, it is necessary to constrain the effective interactions either 
by microscopic many-body calculations for the neutron matter data or the data of nuclei 
with extreme isospin. 

For neutron star matter, the density-dependent interactions DD-MEl and TW-99 favor 
large neutron populations due to their weak p-meson field at high densities. The OV mass 
limits calculated from different EOS are 2.02 ~ 2.81Mq, and the corresponding radii are 
10.78 ~ 13.27 km. The stronger interaction gives a stiffer EOS and a larger mass limit. 
TW-99 and DD-MEl give softer EOS and smaller OV mass limits than those of the effective 
interactions NLl, NL3 and NLSH. The results of TMl are similar to TW-99. After the 
inclusion of hyperons, the corresponding values become 1.52 ~ 2.06 Mq and 10.24 ~ 11.38 
km. 
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TABLE I: The Fermi momenta kp (po = kp/Sir'^), vector densities po, scalar densities ps, scalar 
potentials g^Bf^^ vector potentials gujB^O} the energy densities e^, £ui, ^re, and e^, the 

total energy densities e, nucleon masses m and energies per nucleon Eb/A = e/ p — m for different 
interactions in symmetric nuclear matter at saturation density (the units are in MeV- fm~^ except 
otherwise stating). 





TW-99 


DD-MEl 


GL-97 


NLl 


NL2 


NL3 


NLSH 


TMl 


TM2 


A;F(fm-i) 


1.313 


1.310 


1.313 


1.310 


1.292 


1.300 


1.293 


1.291 


1.251 


/9o(fm-3) 


0.153 


0.152 


0.153 


0.152 


0.146 


0.148 


0.146 


0.145 


0.132 


Ps(fm"3) 


0.143 


0.143 


0.148 


0.142 


0.139 


0.140 


0.138 


0.138 


0.125 


5^BC7(MeV) 


-417.9 


-396.2 


-206.4 


-400.7 


-309.7 


-380.3 


-378.24 


-342.94 


-402.06 


-^aBCrps 


29.871 


28.279 


15.274 


28.450 


21.524 


26.621 


26.099 


23.663 


25.129 




29.871 


28.279 


13.830 


28.187 


20.538 


26.088 


25.544 


22.387 


24.027 




338.7 


316.7 


145.57 


325.4 


242.9 


308.0 


306.6 


274.5 


331.8 


— ^9ojBi^0P0 


-25.910 


-24.069 


-11.140 


-24.689 


-17.688 


-22.832 


-22.378 


-19.901 


-21.897 


£uj 


-25.910 


-24.069 


-11.140 


-24.689 


-17.688 


-22.832 


-22.378 


-19.408 


-21.221 




-0.161 


-0.364 























Eat 


51.981 


48.502 


22.281 


49.378 


35.375 


45.664 


44.757 


39.856 


43.893 




85.401 


87.912 


116.199 


87.207 


96.244 


87.643 


86.786 


91.016 


75.261 




137.382 


136.414 


138.470 


136.585 


131.619 


133.307 


131.543 


130.872 


119.155 


£ 


141.182 


140.260 


141.169 


140.083 


134.470 


136.563 


134.708 


133.649 


120.678 


m(MeV) 


939 


939 


939 


938 


938 


939 


939 


938 


938 


EB/A{MeV) 


-16.25 


-16.23 


-16.32 


-16.42 


-17.02 


-16.25 


-16.35 


-16.265 


-16.16 
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TABLE II: Nuclear matter saturation properties for different effective interactions, including the 
Fermi momenta kp, saturation densities pQ, energies per nucleon Eb/A, effective masses m* and 
m*/m, incompressibility K and symmetric energy coefficients asym- 





A>(rm-i) 






iu*(MoV) 


in* / in 


A'(MoV) 


a,sy,„(M(;V) 


DD-MEl 


1.310 


0.152 


-16.23 


542.7899 


0.578 


244.50 


33.06 


TW-99 


1.313 


0.153 


-16.25 


521.0724 


0.555 


240.00 


32.77 


NLl 


1.310 


0.152 


-16.42 


537.2949 


0.573 


211.15 


43.47 


NL2 


1.292 


0.146 


-17.02 


628.3138 


0.670 


399.17 


45.12 


NL3 


1.300 


0.148 


-16.25 


558.6835 


0.595 


271.73 


37.42 


NLSH 


1.293 


0.146 


-16.35 


560.7559 


0.598 


355.34 


36.12 


TMl 


1.291 


0.145 


-16.265 


595.0626 


0.634 


281.17 


36.89 


TM2 


1.251 


0.132 


-16.16 


535.9376 


0.571 


343.83 


35.98 


GL-97 


1.313 


0.153 


-16.32 


732.6145 


0.780 


240.00 


32.5 



TABLE III: The central densities, energy densities, and pressures, OV mass limits and correspond- 
ing radii for neutron stars for different effective interactions. For interactions DD-MEl, TW-99, 
TMl and GL-97, the same quantities for neutron stars with hypersons are given in the following 

lines r('spcx-ti\'ely. 





Central density 
(fm-3) 


Central energy density 
(xl0^5g/cm^) 


central pressure 
(xlO^^dyne/cm^) 


OV mass limit 
(xMq) 


Radius 
(km) 


DD-MEl 


0.815 


1.852 


7.886 


2.475 


11.903 




0.980 


2.097 


5.780 


2.061 


11.375 


TW-99 


0.970 


2.126 


7.902 


2.195 


11.209 




1.176 


2.422 


6.318 


1.868 


10.853 


TMl 


0.852 


1.881 


5.293 


2.180 


12.054 




1.016 


2.083 


3.071 


1.517 


11.366 


GL-97 


1.045 


2.347 


7.905 


2.018 


10.779 




1.299 


2.792 


6.188 


1.610 


10.242 


NLl 


0.658 


1.529 


7.140 


2.809 


13.137 


NL3 


0.667 


1.544 


7.055 


2.778 


13.081 


NLSH 


0.649 


1.497 


6.682 


2.803 


13.270 
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FIG. 1: (Color online) The effective interaction strengths for a (top), uj (middle) and p (bottom) in 
symmetric nuclear matter as functions of the nucleon density. The shadowed area corresponds to 
the empirical value of the saturation density in nuclear matter ( Fermi momentum kp = 1.35 ±0.05 
fm~^ or density p = 0.166 =b 0.018 fm~^). The curves are labelled from the top to the bottom at 
p=0.15 fm~^ orderly from left to right. 
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FIG. 2: (Color online) Similar as Fig. ^ but for neutron stars matter. The curves are labelled 
from the top to the bottom at p=0.9 fm^^ orderly from left to right. 
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FIG. 3: (Color online) Vector potentials guiB^o (top), scalar potentials gaB(^ (middle) and the sum 
of both (bottom) in symmetric nuclear matter for different effective interactions (as marked in the 
figure) as functions of the nucleon density p. The shadowed area corresponds to the empirical value 
of the saturation density in nuclear matter. The curves are labelled from the top to the bottom at 
/9=0.15 fm~^ orderly from left to right. 
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FIG. 4: (Color online) Vector potentials QujB'^o-, scalar potentials QaSO', isospin vector potentials 
9pBPo,3 and the sum of them in pure neutron matter for different effective interactions as functions 
of the nucleon density p. The curves are labelled from the top to the bottom at p=0.15 fm~^ 
orderly from left to right. 
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FIG. 5: (Color online) Vector potentials g^B^o-, scalar potentials g'o-scr, isospin vector potentials 
9pBPo,3 and the sum of them in neutron stars with neutrons and protons only, for different effective 
interactions as functions of the baryon density p. The curves are labelled from the top to the 
bottom at p=0A5 fm^^ orderly from left to right. 
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FIG. 6: (Color online) Energies per nucleon, Eb /A = e/p—m, in pure neutron matter (upper panel) 
and symmetric nuclear matter (lower panel) for difference effective interactions as functions of 
nucleon density p. The dots in the lower panel correspond to the saturation densities in symmetric 
nuclear matter and the shadowed area corresponds to the empirical value of saturation densities 
in symmetric nuclear matter (density p = 0.166 it 0.018 fm~^ and energy per particle e/p = 
—16.0 lb l.OMeV). The curves are labelled from the top to the bottom at p=0.2 fm~^ orderly from 
left to right. 
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FIG. 7: (Color online) Similar as Fig. |3 but for neutron stars. 
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FIG. 8: (Color online) Particle densities (from top to bottom are respectively n, p, e~ and /i~ 
) in neutron stars for different effective interactions as a function of the baryon density p. The 
corresponding inserts are the same figures in a logarithmic scale at low density (0.05 < p ire? < 0.2). 
The curves are labelled from the top to the bottom at p=0A5 fm~^ orderly from left to right. 
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FIG. 9: (Color online) EOS of neutron stars for different effective interactions (as marked in figure). 
The solid lines represent these with neutrons and protons only, and the dashed lines represent the 
corresponding ones with hyperons included respectively. 




FIG. 10: (Color online) The masses versus the central densities (left panel) and radii (right panel) 
in neutron stars for different effective interactions (as marked in figure). The solid lines represent 
these with neutrons and protons only, and the dashed lines represent the corresponding ones with 
hyperons included. 
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